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Products of positive semi-definite matrices 
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Dedicated to Professor Rajendra Bhatia. 

Abstract 

It is known that every complex square matrix with nonnegative determinant is the product 
of positive semi-definite matrices. There are characterizations of matrices that require two or 
five positive semi-definite matrices in the product. However, the characterizations of matrices 
that require three or four positive semi-definite matrices in the product are lacking. In this 
paper, we give a complete characterization of these two types of matrices. With these results, 
we give an algorithm to determine whether a square matrix can be expressed as the product of 
k positive semi-definite matrices but not fewer, for fc = 1, 2, 3,4, 5. 
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1 Introduction 


Let Mn be the set of n x n complex matrices. In [3], the author showed that a matrix in Mn with 
nonnegative determinant can always be written as the product of five or fewer positive semi-definite 
matrices. This is an extension to the result in [T] asserting that every matrix in with positive 
determinant is the product of five or fewer positive definite matrices. Analogous to the analysis in 
[T], the author of [3] studied those matrices which can be expressed as the product of two, three, or 
four positive semi-definite matrices. In particular, characterization was obtained for the matrices 
that can be expressed as the product of two semi-definite matrices; also, it was shown that any 
matrices not of the form zl, where is not a nonnegative number, could be written as the product 
of four positive semi-definite matrices. Moreover, it was proved in [3l Theorem 3.3] that if one 

\Ti 


applies a unitary similarity to change the square matrix to the form T = 


0 


* 

T2 


G M„ such that 


Ti is invertible and T 2 is nilpotent, and if Ti is the product of three positive definite matrices, then 
T can be expressed as the product of three positive semi-definite matrices as well. It was suspected 
that the converse of this statement is also true. However, the following example shows that it is 
not the case. 
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Example 1.1 Let T = 


Then T = 


■-9 -9 
0 0 

positive semi-definite matrices. However, Ti 
matrices because det(Ti) < 0. 


■9 3 


■ 13 -15' 


T T 

3 2 


-15 18 


1 1 


is the product of three 
[—9] is not the product of three positive definite 


Of course, one may impose the obvious necessary condition that det(Ti) > 0, and ask whether 
the conjecture is valid with this additional assumption. Nevertheless, it is easy to modify Example 
1 by considering 


T(g)l2 


—9/2 —9/2 

O2 O2 


{Ai 0 12)(^2 ® 1^2)(^3 ® -^3) 


using the factorization T = A 1 A 2 A 3 . In the modified example, we have Ti = — 9/2 and T 2 = O 2 . 
By m Theorem 4], Ti = — 9/2 is the product of no fewer than five positive definite matrices. 

In the next section, we will give a complete characterization for those matrices that can be 
written as the product of three positive semi-definite matrices and not fewer. Also, we add an 
easy to check necessary and sufficient condition for invertible matrices that can be written as the 
product of three positive definite matrices. With these results, one can use the Jordan form of a 
given matrix and its numerical range to decide whether it can be expressed as the product of two, 
three, four or five positive semi-definite matrices. 


2 Product of three positive semi-definite matrices 


We will prove the following. 
Theorem 2.1 Suppose 


Ti R 
0 T 2 


such that Ti is invertible and T 2 is nilpotent. 


( 1 ) 


Then T is a product of three positive semi-definite matrices if and only if one of the following holds. 

(a) R 0 or T 2 7 ^ 0. 

(b) R = 0, T 2 = 0, and Ti is the product of three positive definite matrices. 

We establish some lemmas to prove Theorem 12.11 The first one is covered by [TJ Theorem 1]. 

Lemma 2.2 Let A,S€z such that S is invertible. Then A is a product of an odd number of 
positive semi-definite matrices if and only if S*AS is. 

In the next lemma, we need the concept of the numerical range of a matrix A € defined by 

W{A) = {x*Ax : X € C^,x*x = 1}. 

The numerical range is a useful tool in studying matrices. One may see [SJ Chapter 1] for the basic 
properties of the numerical range. 
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Lemma 2.3 Suppose T = 


Ti R 
0 0 , 


pj 


such that Ti G Mm is invertible and R is nonzero. Then T is 


a product of three positive semi-definite matrices. 


Proof. First, we show that there is a px m matrix S such that Ti + RS is the product of three 
positive definite matrices. 

If m = 1, there is 5 € such that Ti + RS > 0. Suppose m > 1 and R has a singular value 
decomposition R = sixiy* + ■ • • + SkXkUl, where si,. ■ ■ ,Sk are nonzero singular values of R, and 
xi,... G C”* and yi,... ,yk G are their corresponding right and left unit singular vectors 
accordingly. Let {ei,..., Cm} be the standard basis for C”^. Take a unitary U such that Uxi = ei. 
Since Ti is invertible, so is Ti = UTiU*. Suppose is the first row of Ti. Let v = ti + ee 2 with 
e > 0 and ri(e) be the m x m matrix obtained from Ti by replacing its first row with v*. Take a 
sufficiently small e > 0 such that v is not a multiple of ei and the matrix ri(e) is still invertible. 
Set S = sf^re^^yiv*U with r > 0 and 6 G [0,27r). Then 

U{Ti + RS)U* = UTiU* + URSU* = fi + re^^eiv*. 

Since v is not a multiple of ei, the rank one matrix eiv* is not normal and so W(eiv*) is an elliptical 
disk with foci 0 and vi with length of minor axis y^||u|P — |uip > 0, where vi is the first entry of 
v\ for example, see O Theorem 1.3.6]. By the fact that the map X i—?■ W{X) is continuous, for a 
sufficiently large r > 0, W (^iv* + e“*®ri/r^ still contains 0 as an interior point for any 9 G [0, 27r). 
Then so does 

W (fi + = re*® W (^eiv* + • 

In addition, the value r can be chosen so that r > | det(Ti)/det(T'i(e))|. Now by the linearity of 
determinant with respect to the first row, 

det(ri + RS) = det(C/(ri + RS)U*) = det(fi + re*®eiu*) = det(fi) + re*® det(fi(e)). 

Since |det(Ti)| < |re*® det(Ti(e))|, there is 0 G [0, 27r) such that det(Ti + RS) > 0. By [H Theorem 
3] (see also Proposition 12.4p . Ti + RS is a product of three positive definite matrices. 

Finally, note that 


'Im S* 

'Ti R 

Im 

o' 


'Ti + RS R 

0 Ip 

0 Op 

_s 

Ip 


0 Op 


By [21 Theorem 3.3], T is a product of three positive semi-definite matrices, and so is T by Lemma 

1121 □ 

Proof of Theorem \2.1\ . Suppose T is the product of three positive semi-definite matrices. If 
R 0 or T 2 0, then we are done. Else, T = Ti © Op is the product of three positive semi-definite 
matrices. By [3l Proposition 3.5], Ti is the product of three positive definite matrices. 
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To prove the converse, we consider the following three cases. 

Case 1. Suppose R ^ 0. We use induction on p, the size of T 2 . If p = 1, the result follows from 
Lemma f2.3l as T 2 = [0]. Assume the result holds for T 2 with size at most p—l- Since T 2 is nilpotent, 

[221 T22 


without loss of generality, we may assume that T 2 is upper triangular with T 2 = 


0 0 


, where 


T 21 G Mp_i. Write R = [i?i R 2 ] where is m x (p — 1). If Ri / 0, then by induction, the 
matrix 

U J 21 

Theorem 2.2] (see also Proposition 13.1|) . we may assume that both Pi and P 2 are invertible. Let 
'R 2 


X = 


T 22 


is a product of three positive semi-definite matrices, say, P 1 P 2 P 3 . Further, by [3l 
3 also Proposition ED, we may assume t 
with size (m + p — 1) x 1. Then for any e > 0, 


'Pi o' 

P 2 

ePf^A' 

1 

0 

£ 


'P 1 P 2 P 3 X 

0 0 

e(Pf^A)* 

1 

0 e-i 


0 0 


Ti Ri R2 

0 r2i T 22 

000 


= T. 


Clearly, Qi = Pi® [0] and Qs = P 3 ® [e are positive semi-definite matrices. Now one can choose 


a sufficiently small e > 0 so that Q 2 = 


P 2 

e(Pf^A)* 1 


is also positive semi-definite. Thus, 


T is a product of three positive semi-definite matrices QiQ 2 Q 3 - 

Now suppose Ri = 0. Then the (m -|- l)th column of T is a zero column. By interchanging the 
{m + l)th and the last indices, one can see that T is permutationally similar to 


Ti R O' 
0 221 0 

0 T 22 0 


where R is m x (p — l), T 21 is (p — 1 ) x (p — 1 ), and T 22 is 1 x (p — 1 ). 


Notice also that R is nonzero and T 21 is nilpotent. By induction, 


is a product of three 


Ti R 

0 T 21 

positive semi-definite matrices P 1 P 2 P 3 . By [3l Theorem 2.2], we can further assume that both P 2 
and P 3 are invertible. Let T = [O 2 ^ 22 ] with size 1 x (m + p — 1). Then for any e > 0, 


Ti R O' 
0 f2l 0 
0 T 22 0 


'Pi 0 ■ 

P 2 

<YP^Y 

P 3 0 


P 1 P 2 P 3 0 

0 e"^ 

.<YP3^) 

1 

0 0 


Y 0 


Again, one can choose a sufficiently small e > 0 such that all three matrices in the left side of 
the above equation are positive semi-definite. Thus, T is permutationally similar to a product of 
three positive semi-definite matrices, and hence T can also be written as a product of three positive 
semi-definite matrices. 

Case 2. Suppose R = 0 and T 2 is nonzero. Without loss of generality, we may assume that T 2 
is upper triangular with zero diagonal entries while the first row of T 2 is nonzero. Let Z be the 
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p X m matrix with 1 at the (1, l)-entry and zero elsewise. Then Z*T 2 ^ 0 and T 2 Z = 0. Let 
r/m O' 


5 = 




Then 


'Im Z* 

Tl 

O' 

Im 

o' 


'Ti + Z*T2Z 

Z*T2 


Ti 

Z*T2 

. 0 Ip. 

0 

T2_ 

z 

Ip. 


T 2 Z 

T 2 _ 


0 

T 2 _ 


Since Z*T 2 7 ^ 0, the result follows from Case 1 and Lemma 12.21 


Case 3. Suppose i? = 0 and T 2 = 0. If Ti is the product of three positive definite matrices, then 
T = Ti 0 Op is the product of three positive semi-definite matrices. □ 

Note that Theorem 12.11 depends on checking an invertible matrix is the product of three positive 
definite matrices. Such conditions are given in [H Theorem 3]. We restated the result in terms of 
the numerical range in the following proposition, which is based on the discussion in [H Theorem 
3 and Fact 3.2]. 


Proposition 2.4 Let T € be such that det(T) > 0. Then T is the product of three positive 
definite matrices if and only if one of the following holds. 

(a) W{T) contains 0 as an interior point. 

(b) W{T) contains a positive number, and the arguments of the eigenvalues ofT can be arranged 

as: —vr < 61 < ■ ■ ■ < On < such that — 0 - 

Note that in [H p. 88 ], the author required in condition (3.6b), corresponding to our condition 
(b), that all real eigenvalues of T are positive, which is ensured by our assumption that 9j G (—vr, vr) 
for all j. 


3 Determining the number of factors 

In this section, we describe an algorithm to determine the smallest number of positive semi-definite 
matrices whose product equals a given A G with nonnegative determinant. 

We first present the following theorem providing some easy tests for a matrix A to be the 
product of two positive semi-definite matrices. The equivalence of conditions (a), (b), (c) were 
given in O Theorem 2.2]. We include a short proof, which is different from that of Wu. 

Proposition 3.1 Let A be a square matrix. The following are equivalent. 

(a) A is the product of two positive semi-definite matrices. 

(b) A = BC, where B, C are positive semi-definite matrices such that B or C can be assumed to 
be invertible. 

(c) A is similar to a nonnegative diagonal matrix. 
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(d) A is unitarily similar to an upper block triangular matrix such that the diagonal blocks are 
scalar matrices corresponding to distinct scalars. 

(e) The minimal polynomial of A only has simple nonnegative zeros. 


Proof. The equivalence of (c), (d), (e) are clear. 

(c) ^ (b): A = S-^DS = S-\S-^y{S*DS) = S-^D{S-^)*{S*S). 

(b) (a): Trivial. 


be 


(a) ^ (c): Suppose A = BC, where B and C are nxn positive semi-definite. Let U be unitary 
such that U*BU = Bq (B 0^, where Bq G M^-k is positive definite. Let U*CU = 

[C21 C22 

such that C 22 £ ^k- Assume V is unitary such that V*CiiV = Cq © 0^ for a positive definite 
matrix Cq G Mn-k-e- We may replace U by U(V (B I) and assume that Cn = Cq © 0 ^. Because C 

'Co 0 Cl 


is positive semi-definite, U*CU = 


U*AU = 


0 Of 0 


Now 


1 


Bo^ 0 


0 4 


1 


Bo^ 0 


0 

- 


[ 

Cj* 

0 

0 


0 

0 

- 

0 Cq-^Ci] 

I 0 

0 4 



. Then 


Co 0 Cl 
0 Of 0 
Cl* 0 C 22 J 


Bo 0 
0 Ofc 


Co 0 Cl' 
0 Of 0 

0 0 Ofc 


Bo 0 
0 Of; 


Co 0 Cl' 
0 Of 0 
0 0 Of; 


T 0 

0 I 
0 0 


-Co-'Ci' 
0 
4 




0 


Lo 

4_ 


Bo 0 
0 Of; 



0 

0 


1 

1 ^ . 



0 Of 0 


Bf 0 

cq 

II 

Co 0 

- 

_ 0 0 0f;_ 


0 4 

V 

0 Of 


Bf 1 © Of;. 


Therefore, A is similar to 


Co 0 
0 Of 

similar to a nonnegative diagonal matrix. 


R2 


R2 


©Of;, which is positive semi-definite and is (unitarily) 

□ 


Now, we are ready to present an algorithm to check whether a matrix A G with det(A) > 0 
can be written as a product of k positive semi-definite matrices, but not fewer, for k = 1,2, 3,4, 5. 


Algorithm 3.2 Let A G A4 with det(yl) > 0. 

If A = aln such that a ^ [0,oo), then A can be expressed as a product of five positive semi- 
definite matrices, but not fewer. Otherwise, apply a unitary similarity to A to get an upper 
triangular matrix 


Ti R 
0 T 2 


such that Ti is invertible and T 2 is nilpotent. 
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( 1 ) IfT is a nonnegative diagonal matrix, then A is itself a positive semi-definite matrix. 

(2) Condition (1) does not hold, and A satisfies any one of the equivalent conditions in Propo¬ 
sition ixn Then A can be expressed as a product of two positive semi-definite matrices, but 
not fewer. 

(3) Suppose (1) and (2) do not hold. Then A can be expressed as a product of three positives 
semi-definite matrices, but not fewer, if any of the following holds. 

(3.a) i? or r 2 is nonzero. 

(3.b) Both R = 0 and T 2 = 0, and Ti is the product of three positive definite matrices. 

In (3.b), the invertible matrix Ti is a product of three positive definite matrices if one of the 
following holds. 

(i) det(ri) > 0 and IT(ri) contains 0 as an interior point, 

(ii) 0 is not in the interior of W{Ti) and W{Ti) contains a positive number and ~ 
where —vr < 0i < ■ ■ ■ < 9^ < tt are the arguments of the eigenvalues of Ti. 

(4) Suppose conditions (1), (2), (3) do not hold, i.e., T = Ti © Op such that neither (i) nor 
(ii) holds for the upper triangular matrix Ti. Then A can be expressed as a product of four 
positive semi-definite matrices, but not fewer. 
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